UTHEP-296 
February 1995 



Deconfining Chiral Transition in QCD 

on the Lattice* 



ON 
D 

P-H 



K. Kanaya^ 

Institute of Physics, University of Tsukuba, Ibaraki 305, Japan 



Abstract 



{Nl The deconfining chiral transition in finite-temperature QCD is studied on 

the lattice using Wilson quarks. After discussing the nature of chiral limit 
with Wilson quarks, we first study the case of two degenerate quarks Np = 2, 
q and find that the transition is smooth in the chiral limit on both N t = 4 and 6 

CN lattices. For Np = 3, on the other hand, clear two state signals are observed 

^ for m q < 140 MeV on N t = 4 lattices. For a more realistic case of Np = 2 + 1, 

ON i.e. two degenerate u and d-quarks and a heavier s-quark, we study the cases 

m s ~ 150 and 400 MeV with m u = ~ 0: In contrast to a previous result 
with staggered quarks, clear two state signals are observed for both cases, 
Qh suggesting a first order QCD phase transition in the real world. 
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1 INTRODUCTION 



The study of QCD on the lattice has shown that QCD has a finite temper- 
ature phase transition/crossover at T ~ 100 - 200 MeV. In the light quark 
limit (chiral limit), the chiral symmetry, that is broken spontaneously at low 
temperatures, is resored when T gets sufficiently high. In the heavy quark 
limit (pure gauge limit), existence of a deconfining phase transition separat- 
ing the low temperature confining phase and the high temperature deconfin- 
ing phase is established. Studies on the lattice for intermediate quark masses 
indicate that the chiral transition in the light quark limit and the deconfin- 
ing transition in the heavy quark limit are two ends of a single transition 
line with a probable crossover part in the middle. This transition/crossover 
plays a primary role in the studies of the quark-gluon- plasma state in heavy 
ion collisions and in the evolution of early Universe. A decisive information 
for these studies is the order of the transition. Therefore it is important to 
pin down the nature of the deconfining chiral transition with the values of 
parameters corresponding to the real world. 

Previous investigations based on universality arguments as well as numer- 
ical studies on the lattice show that the order of the transition in the chiral 
limit depends on the number of flavors Np: With degenerate Np quarks, the 
chiral transition is first order for Np > 3, while a continuous transition is 
suggested for Np = 2. Because the s-quark mass m s is of the same order 
of magnitude as the transition temperature, the order of the transition may 
depend on the value of m s sensitively. We need to include the s-quark to 
make a prediction for the real world. 

In this report, we study the finite temperature transition on the lattice 
using Wilson fermions [1] as lattice quarks. The action is given by 





S g = -£ E Tr(P p + Pj) 



(2) 



plaquettes 

N F 



(3) 



f=l x,y 
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where j3 = Q/g 2 with gauge coupling g, 

4 

~ K f ~ lnWx^x + fx, y + (1 + j^Uy^Sx, y + ^}, 

and is the hopping parameter to tune the quark mass for the f-th flavor. 
With SU(3) gauge group lattice spacing a is a decreasing function of p and 
the temperature is given by T = l/(N t a(p)) with N t the lattice size in the 
temporal direction. 

In the next section, we first discuss the chiral limit with Wilson quarks. 
We then present the results of our simulations for the case of degenerate 
quarks with Np = 2 and 3. Finally we discuss a more realistic case of 
Np = 2 + 1: two degenerate light u and d-quarks and a heavier s-quark. 
Comparison with previous results using staggered quarks is also made. We 
mainly perform simulations on 8 2 x 10 x N t and 12 3 x N t lattices with N t = 
4. For Np = 2 we also study the case of N t = 6. Preliminary reports of the 
resirits presented here can be found in refs. [2, 3, 4, 5] and [6]. More detailed 
reports will be available soon. 

2 Chiral limit 

Before studying the finite temperature transition, let us discuss the definition 
of chiral limit with Wilson quarks. In the continuum limit p = oo, it is 
easy to see that the quark masses vanish at K = 1/8 = Kq(P = oo) as 
m q oc 1/K — 1/Kc- When the system is in the chirally broken phase, we 
expect ml oc m q and Kc can also be defined as the point where m\ vanishes. 

2 

Let us denote Kc defined through m q as K c " and that through m\ as K™* 
where a distinction among these two definitions is needed. For finite lattice 
spacings a > 0, however, Wilson's action (3) breaks the chiral symmetry 
explicitly through the Wilson term and the notion of chirality as well as the 
definition of Kc becomes unclear. 

Nevertheless numerical simulations at finite p's show that, in the low 
temperature confining phase (with iV^ < 6), vanishes linearly in 1/K at 
finite Kc(P)- 1 Accordingly, a careful perturbative study of Wilson quarks at 
large /?'s show that, when renormalization counter terms are chosen properly, 

1 In the space of coupling parameters, it is difficult to simulate the points where is 
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axial Ward identities hold also on the lattice to the price of 0(a) correction 
terms [7]. We can therefore consider, for example, a relation: 

2m q < | P 1 7r >= -m % < | A 4 | tt > +0(a) (A) 

to define m q where P is the pseudoscalar density and A 4 the fourth com- 
ponent of the local axial vector current. Here, the chiral limit Kq can be 
defined where m q vanishes. Kq gets a /3-dependence through counter terms. 

Independently Tsukuba group [8] proposed to define quark mass through 
an axial Ward identity neglecting 0(a) terms in (4). Maiani and Martinelli 
[9] also studied m q defined similarly using axial Ward identity. Numerical 
studies of this m q for intermediate (5 ((3 = 5.85 with quenched QCD [10] and 
P = 5.5 with Np = 2 full QCD [11]) show that (i) m q vanishes linearly in 
1/K at finite K = K™", (ii) the value of m q (therefore the value of K™" also) 
is almost independent of T and the phase, and (hi) K™" is consistent with 

2 

K^"" determined with m 2 in the confining phase (Fig. 1). Therefore we have 

if™' (high T) ^ K™ q (\owT) ~ if™* (low T). 

These nice features of our m q suggest that 0(a) chiral violation are small in 
(4) down to these values of j3. We also note that we can simulate the system 
in the deconfining phase even just on the i^c-line without encountering the 
critical slowing down. Our m q therefore provides us a numerically easy way 
to determine Kq for f3 > 5.5. 

2 

Equivalence of K™" and if™ x can be shown also in the strong coupling 
limit (3 = [12]: In his pioneering paper [1], K. Wilson carried out a strong 
coupling calculation of meson masses for T = without quark loops. We can 
extend this calculation to our m q and we get 

(1 - 16K 2 )(1 -4K 2 ) 
cosh(m v a) = 1 + — — - (5) 

AK 2 sinh(m T a) 
2m " a = m * a i-AK 2 co^(m„a) (6) 

Our numerical data for Np = 2 at (3 = agrees with these formulas well 
as shown in Fig. 2. (Results for Np = 18 is presented in ref. [12] where the 



very small due to the critical slowing down. Here, -K - ™* is defined by a linear extrapolation 
of m\ obtained in the confining phase on a lattice as large as possible. 
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agreement in this case is shown. The rho meson mass, the nucleon mass, and 
the delta mass also agree with corresponding mass formulas.) We note that 
both m v and m q vanish at the same K = 0.25 and oc m q near this Kc- 

For < /3<5.3, however, an additional complexity comes on the stage: 
m q shows a strange dip when we increase K over the deconfining chiral 
transition/crossover point and m q in the deconfining phase is no more 
independent of T for these /3's [2, 12, 13, 14]. This strange behavior can be 
attributed to an effect of 0(a) chiral violation of Wilson fermions. In the 
confining phase, on the other hand, no such strange behavior is found. This 
interpretation including the smallness of 0(a) effects in the confining phase 
is confirmed by our recent study with an RG improved action [4, 5]. For 
these /3's, we define Kc by a linear extrapolation of or m q obtained in 
the confining phase on a lattice as large as possible. 

Fig. 3 summarizes the numerical results of Kc for Np = 2 including the 
data from other collaborations [13, 14, 15, 16]. The sharp change of Kc 
around p « 5.0 signals the crossover from the strong coupling region to the 
weak coupling region. Our data for p < 4.3 are obtained mainly on an N t = 4 

2 

lattice. Small difference in the value of K™* and K c " there (at most of the 
order of 0.01) will be a result of errors in the extrapolation to = and 
2m q = 0, and also of the O(a) chiral violation by Wilson fermions — our data 
is not accurate enough yet to decide the main cause among these possibilities. 

2 

Taking the errors due to the difference of K™* and K™ q into account, we 
find that N t (as well as Np) dependence of Kc is small for p < 4.5. 

In summary, in spite of the explicit chiral violation, Wilson quarks show 
rather clean chiral behaviors on the lattice. Although some care is required 
in handling numerical results of m q etc. in the deconfining phase at small 
/3's, most effects of chiral violation seem to be absorbed by renormalizations 
such as a shift of Kq and Z-factors for axial currents etc. 

3 N F = 2 

The location of can be easily identified numerically by rapid changes of 
physical observables such as plaquette, Polyakov loop, hadron masses, and 
m q . Close to the chiral limit Kc, however, the critical slowing down in the 
confining phase makes the simulation practically impossible. In numerical 
simulations of QCD, this critical slowing down appears as an rapid increase 
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of the number 7V; nv of CG iterations to invert the quark matrix, closely con- 
nected with the appearance of very small eigen values of quark matrix [11]. 
As discussed before, on the other hand, A^; nv remains small in the deconfining 
phase even just on the i^c-rme. We can therefore study the chiral transition 
by simulations along the iQ^-lme in the deconfining phase: Decreasing (3 from 
the deconfining phase, the location of the chiral transition, i.e. the crossing 
point [3ct of Kt an d Kq lines, can be determined by a sudden increase of 
iVi nv . The value of (3ct determined by this method is consistent with an ex- 
trapolation of the i^-line determined by the conventional method as shown 
in the phase diagrams given below. The nature of the transition there can 
also be studied by monitoring physical observables on the i^-line. 

In Refs. [3, 6, 17] this method is applied to the case of Np = 2. Phase 
diagram for Np = 2 is shown in Fig. 4. Time history of A^; nv on the i^c-line 
is given in Fig. 5 for N t = 4. Time histories of plaquette, Polyakov loop, m v 
etc. at the same points confirm that the system is developing into a confined 
state for j3 < 3.9. Therefore we conclude that j3 C T - 3.9 - 4.0 for N t = 4. 
Similar study on an N t = 6 lattice gives (3ct — 4.0 - 4.2 and a study on a 
18 3 x 24 lattice gives per ~ 4.5 - 5.0. 

Our result of m 2 on the i^c-line is shown in Fig. 6. For both N t = 4 and 
6, we find that, when we decrease (3 toward Act, decreases rapidly and is 
consistent with zero at (3ct- This suggests that the transition is smooth in 
the chiral limit, in accord with the residt of an universality argument [18] and 
finite-size scaling studies with staggered quarks [19, 20] suggesting a second 
order transition. Off the chiral limit, m q > 0, we may therefore expect that 
Kt is crossover. 2 

2 Concerning the nature of the transition on the ify-line at intermediate values of 
K, MILC collaboration reported recently an unexpected phenomenon that the transi- 
tion/crossover becomes once very strong when we increase m q , while it becomes weak 
again when we further increase m q [13, 14]. For N t = 4, Kt is very strong at K ~ 0.18. 
For 7V t = 6 they even found first order signals at several iT's. Consulting the phase di- 
agram Fig. 4, we note that these values of K just correspond to the region where the 
i^c-line shows a sharp bend due to the crossover phenomenon between weak and strong 
coupling regions. The ii^-line that deviates from the i^c-line at (3ct gets very close to the 
i^c-hne again there. It is therefore plausible that the strong Kt is a lattice artifact in the 
crossover region for these 7V t caused by this unusual relation of Kq and Kt lines. Close 
to the continuum limit with a large enough N t , the _K"c-hne is smooth and the ify-line 
simply gets apart from it when we increase m q . 
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4 Np = 3 



We apply the same method to study the chiral transition for Np = 3 and we 
find a two-state signal at (3ct — 3.0 at N t = 4 [3, 17]: When we start from a 
deconfining configuration prepared with a larger f3, it remains stable in the 
deconfining phase with small A^ nv . However, when we start from a mixed 
initial configuration, 7V; nv increases rapidly with MC-time and simultaneously 
plaquette, Polyakov loop, and decrease rapidly, showing that the system 
is developing into a confining state. This implies a first order chiral transition 
for Np = 3. Similar first-order signal is obtained also for Np = 6 [2]. These 
results are consistent with the prediction based on universality [18]. Our 
results of on the i^-line is summarized in Fig. 7. 

We now extend the study to finite m 3 's by increasing (3 on the i^-line. 
Our phase diagram for Np = 3 is given in Fig. 8. We find clear two-state 
signals for (3 = 4.0, 4.5 and 4.7 on both 8 2 x 10 x 4 and 12 3 x 4 lattices 
(Fig. 9). However, when we further increase f3 > 5.0, we find no clear signs 
of met ast ability. The critical value of the quark mass m™*, up to which we 
get a clear first order signal on these lattices, is bounded from below by m q 
for = 4.7. We find mf x a > 0.175(2) and (m T /m p ) crit > 0.873(6). 

In order to convert the critical quark mass into physical units, we use 
the results of m p a and m q a for these /3's. When we plot the masses in the 
confining phase as a function of 1/K — 1/Kc, they are almost independent 
of p, Np and N t for /3<4.7. Fig. 10 shows our results for Np = 2 and 
3. Identifying m p (K c ) with 770 MeV we get a ~ 0.8 GeV -1 for /3 < 4.7. 
It should be noted that the vector meson mass "m p " for m q ~ 150 MeV 
is consistent with = 1020 MeV ((f) ~ ss due to the approximate "ideal 
mixing"), indicating that our definition of m q is close to that of the current 
quark mass. 

From this value of a we get m" 1 * > 140 MeV. We note that our m" 1 * is 
much larger than that obtained previously with staggered quarks at N t = 4 
with spatial 8 3 - 16 3 lattices [21, 22]: two state signals are observed for m q a = 
0.025, while no clear met ast abilities are found for m q a = 0.075. Using a result 
for meson masses at corresponding /3 [23] (the values for Np = 4, because they 
are only available) we obtain m™* ~ 12 - 38 MeV and (m T /m p ) crit ~ 0.42 - 
0.58 for staggered quarks. 
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5 N F = 2+1 



We now study the case of Np = 2 + 1: m u = < m a . If m a 3> T c then s- 
quark will decouple from the dynamics relevant to the deconfining transition 
and we expect a smooth deconfining transition similar to the case of Np = 2. 
On the other hand, when m s <C T c , the results for Np = 3 implies that 
the transition is first-order. Realistic value for m s seems to be of the same 
order as T c . Therefore we should include the dynamics of s-quark properly. 
In this connection, we note that, when m s < m™* (the critical quark mass 
for the case of degenerate Np = 3), we can expect a first-order deconfining 
transition, while when m u ^ > m™*, we should expect a crossover. Recall 
the results for m q mt discussed in the previous section: m q mt > 140 MeV with 
Wilson quarks and m™* ~ 12 - 38 MeV with staggered quarks. Our large 
m q ni opens the possibility of a first order deconfining transition at quark 
masses corresponding to the real world. 

In order to see what happens with realistic quark masses, we study the 
cases m s ~ 150 MeV and 400 MeV. Corresponding values of K s are obtained 
from Fig. 10. For u and d-quarks we fix m u = ~ 0. We apply the method 
which we used in previous sections to study the chiral transition: Keeping 
K u = Kd = Kc and K s at the point corresponding to the case of m s ~ 150 
or 400 MeV, we decrease j3 until we hit Kt and study the nature of the 
transition there. Fig. 11 is an example of time history we get. We find two 
state signals on an 8 2 x 10 x 4 lattice for both m s = 150 and 400 MeV and 
also on a 12 3 x 4 lattice for m s = 400 MeV. Our results of on the i^c-line 
are shown in Fig. 7. 

With staggered quarks, Columbia group [22] reported that no transition 
occurs at m u a = m&a = 0.025 and m s a = 0.1 (m u ^ ~ 12 MeV, m s ~ 50 
MeV using the value of a from [23]). Although a direct comparison with 
our result is difficult because of the difference in simulation parameters, this 
implies a discrepancy among lattice quark formulations for these /3's and/or 
a sensitive dependence of the nature of on the value of light quark masses 
m u and m^. 
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6 CONCLUSION 



We have studied the nature of the finite temperature transition with Wilson 
quarks for Np = 2, 3 and 2+1. For Np = 2 the chiral transition is smooth on 
an N t = 6 lattice in accord with our previous result at N t = 4. For Np = 3 
at N t = 4, clear two state signals are observed on the Kf-lme for m q < 140 
MeV. For N F = 2 + 1 we have studied the cases m s ~ 150 and 400 MeV 
with m u = rrid — ? an d we have found two state signals for both cases. This 
suggests a first order finite temperature transition in the real world. 

I am grateful to my collaborators Y. Iwasaki, S. Kaya, S. Sakai, and 
T. Yoshie, with whom the studies reported here are done. Simulations are 
performed with HITAC S820/80 at KEK and with QCDPAX and VPP- 
500/30 at the University of Tsukuba. I thank members of KEK and the 
other members of QCDPAX collaboration for their support. This work is in 
part supported by the Grant-in- Aid of Ministry of Education, Science and 
Culture (N0.O6NPO6OI). 
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Figure 1: m 2 a 2 and 2m q a for JVp = 2 and /3 = 5.5 on an 8 2 x 20 x N t lattice 
with N t = 4 (filled symbols) and 8 (open symbols) [11]. Straight lines are 
the results of a linear fit for 2m q (N t = 4) and m^(N t = 8). For this (3 and 
these values of K, N t = 4 means that the system is in the high temperature 
deconfining phase, while N t = 8 means the low temperature confining phase. 
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Figure 2: The same as Fig. 1 for (3 = on an 8 2 x 10 x 4 lattice. Errors 
are smaller than the size of symbols. Solid curves are the results of a strong 
coupling calculation (6). 



12 




Figure 3: The line of Kc for Np = 2. Filled symbols are for K™* and open 
symbols are for K™ q . Circles are for our data. Our K^ q data at (3=6.0 and 
10.0 are determined in the deconfining phase. 
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Figure 5: Time history of iV inv on the K c -\ine for N F = 2 on an 8 2 x 10 x 4 
lattice. 
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Figure 6: m^a 2 on the i^c-line for N t = 4 and 6 with spatial 8 2 x 10 lattice. 
The locations of the crossing point (3cT(N t ) are shown by shaded bars. 
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Figure 7: m^a 2 on the i^c-line for Np =6, 3 and 2+1. The filled symbols are 
for a 8 2 x 10 x 4 lattice, while the open ones are for a 12 3 lattice. The lower 
points at (3ct represent the upper bounds for m 2 obtained by time-histories 
from mixed initial configurations. 
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Figure 8: Phase diagram for Np = 3 at N t = 4 near Kc- Kr-line for Np = 2 
is also plotted. 
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Figure 9: Time histories of the plaquette for Np = 3 at j3 = 4.7 on a 12 3 X 4 
lattice. 
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Figure 10: m p a and 2m q a in the confining phase. Open symbols are for 
N F = 2, (3 = 3.0, 3.5, 4.0, 4.3, and 4.5 on an 8 2 x 10 x 4 lattice. Filled 
symbols are for N F = 3, (3 = 4.0, 4.5 and 4.7 on 8 2 x 10 x 4 and 12 3 x 4. 
Kc(fl) is determined by m 2 and m q for Np = 2, with errors taking into 
account the difference between two definitions. 
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Figure 11: Time histories of the plaquette for Np = 2 + 1 with m a ~ 400 
MeV on a 12 3 x 4 lattice. 
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